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The theory of class struggle is modelled within the framework of statistical physics. Dichotomous 
spin dynamics on a pyramid-shaped hierarchical structure are examined (akin to the Cayley tree). 
A "head node" is placed at the apex. The system embodies a number of "classes" , corresponding 
to different levels of the hierarchy. A class is comprised of nodes that are equidistant from the 
head. Weighted links exist between nodes from the same and different classes. We study the effect 
of these weights on the dynamics. The spin (hereafter, "state") of the head node is fixed, and it 
imposes its state on the rest of the hierarchy. Necessary conditions so that the classes eventually 
repudiate or acquiesce in the state imposed by the head node are found. The results show that, 
to reach unanimity across the hierarchy, it suffices for the head node to make the bottom-most 
class adopt the same state. Then the rest of the hierarchy will inevitably comply, regardless of the 
inter/intra class link configurations. Hence the role of the bottom-most class is signified. This also 
sheds light on the importance of mass media as a means of synchronization between the top-most 
and bottom-most classes. 



I. INTRODUCTION 

A recent infusion is gathering pace in network science 
whose strands span many sociological and economic in- 
quiries. Network approach has been around for a while 
to capture the structural facets of social networks [THE] ■ 
The approach has been more conceptual than mathemati- 
cal. Extraction of macro properties of the social structure 
from the micro behavior is a linchpin of these efforts. 

The study of micro-macro properties of networks and 
network formations took off in Physics with the seminal 
works [7] and also [5] (tangent to [S] and [TO]). 

Statistical mechanics offers it's own treatment of the 
micro-macro entwinement. The study of opinion dynam- 
ics on graphs is closely knit with the simplest binary 
spin states in statistical mechanics (TTJ [12] (envisaged by 
Lenz [T3] and implemented by Ising [T3], arguably the 
most seminal work of Statistical Mechanics). The Ising 
model has been generalized to various topologies [151119) . 
It was extensively applied to fields outside physics. Here 
we give a few examples. It was used to model strikes on a 
plant j2U| ■ In [23 HH the Ising model was used to model 
binary opinions and social influence. In [23j it was applied 
to empirical election results. The Ising model was applied 
to immigration processes in [23]. In [35] it was used to 
model traders in economic markets. In j^H] it is employed 
to model the fragmentation of Yugoslavia. In [27] it was 
used to model tax evasion dynamics. In [28] it was used 
to model the strength of social ties in a social structure. 
The Ising model is used for graph community detection 
in [29Tl31| . In [32] it was used for image segmentation. 

In [33[ 133] j the Ising model has been endowed with dy- 
namics. The spin at each cite interacts with other spins 



and can change in time, and the time evolution of the sys- 
tem is studied. Another kinetic generalization of the Ising 
model is called the voter model, which is solved on arbi- 
trary graphs in [35]. In [36], |37] three states have been 
envisioned and the dynamics has been solved. The q-voter 
model is introduced in |38j . where a cite adopts the state 
of q of its randomly selected neighbors, given that they 
are unanimous. In |39) . each node besides having a binary 
state, has a continuous number that models the influence 
it receives from neighbors. This number is updated based 
on the behavior of neighbors. In [40], spins are endowed 
with inertia. The longer it has been since their spin hasn't 
changed, the less likely it gets for them to flip. Stubborn 
cites, whose spin never change are envisaged in |41H43j . 
Dynamic spin systems have also been employed to study 
opinion dynamics in groups [33H3Z], marketing [3E1I3HL vot- 
ing models [5UH54"] . etc. 

Generally, the Ising framework helps extract the macro- 
scopic properties of the system under certain prescribed 
interaction schemes in the micro level. This is the feature 
that has ignited efforts to further apply this framework to 
contexts outside Physics. 

In this contribution, nodes form a hierarchical structure. 
A "head node" at level y = 0, connects to q nodes at y = 1, 
which in turn are connected to q other nodes at y = 2, 
and so on (see Figure [T]) . The distance of each node to 
the apex is called its level. We also dub "class" , the set of 
nodes with the same distance to the head. Links also exist 
between nodes of the same class, whose number (fraction) is 
constant in the first (second) model. In the last section, we 
consider arbitrary dependence of the number of intraclass 
links on the cass level. 

In the first part we solve the model on a hierarchy that is 
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FIG. 1: The first four levels of a hierarchy with q = 3. The links whithin each class are not shown in the figure. The number of 
these links are modelled in two different ways in the sections below. First it is assumed to be a constant number, then a constant 

fraction of the possible links. 



semi-infinite (limited from one end) . This provides insight 
on cases where the hierarchy is very long and can be locally 
assumed to be infinite from below. In the second part, a 
finite hierarchy is assumed. In this case, the bottom-most 
class, (which in the Marxian jargon is called "the universal 
class") [SHIES] , is situated at y = L. In all sections but 
the last one, we use a continuous approximation in the y 
domain. But in the last section, we solve for the steady 
state in the general discrete case. 

Each node has a state s = ±1. The head node is fixed at 
s = +1 permanently, and others have aleatory states at the 
outset. We seek conditions necessary for the domination 
of the head node. 

Our findings highlight the importance of the role of the 
universal class. We find that if the universal class is fixed 
aligned with the state of the head node, the middle classes 
will conform irrespective of the density of interclass and 
intraclass links, which merely change the speed of the dif- 
fusion. Alternatively, if the universal class constantly op- 
poses the head, fragmented sets of states emerge in the 
middle. Thus it suffices for the head node to align the 
viewpoint of the universal class as planned, and the rest 
will follow. This highlights the importance of mass media. 

The results resonate with Marx's theory of class strug- 
gle 55-57 . Marxian classification embodies only capi- 
tal and property ownership. This segregates the system 
into haves and have-nots, or the bourgeoise and the pro- 
letariat (and a minority in the middle, i.e., petit bour- 
geoise). This membership leads to communal interaction 
due to a common belief of affiliation, and consequently, 
convergence of actions. This convergence of behavior is 
the rationale behind the assumption made in this paper, 
that nodes from the same class have the same state. This 
convergence shapeshifts into an ever-increasing class in- 
terest and consciousness. It then results in class conflict 
and antagonism. This further intensifies the pressure im- 
pinged upon the proletariat, eventuating in a revolution. 
Weber's classification adds other factors, such as status 
groups, yielding an at least two-dimensional classification, 



altering the consequences [55HST] . 

In this paper, it is tacitly assumed that there is a unique 
classifier to the society, based on which, the levels are as- 
signed. Our analysis is conducted under the mean-field 
assumption, due to the common interests that members 
share within a class. A pyramid-shaped structure is con- 
sidered for the hierarchy as described, and the links that 
tie classes together are assigned values for strength, corre- 
sponding to the weight of influence they carry. 

Conformity to a consensus in the group majority and 
obedience to authority is evident in social interactions [58 - 
R)3"] . The theme employed in this paper for the dynamics 
of nodes' states mimics this behavior. Perhaps the sim- 
plest scenario would be that, at each timestep, each node 
finds the fraction of its neighbors who currently agree and 
oppose with it, and with those probabilities follows either 
of the states. So if, for instance, all adjacent nodes are at 
state +1, the node definitely follows. However if 90% of 
them are at state +1, there is a 10% chance that the node 
will adopt —1. This is the essence of the conventional voter 
models [35], which differs from our problem, where links 
have weights. Each node assigns different weights to those 
with whom it shares the class, and those above and below 
it. Different conditions of the relative magnitude on these 
weights deteriorate or facilitate dominance of the states of 
the head node, as shall be clarified. 



II. FIRST MODEL: CONSTANT NUMBER OF 
INTRA-CLASS LINKS 

We situate the head node at level y = 0, which is con- 
nected to q nodes at level y = 1, and each of those nodes 
are connected to q > 1 new nodes at level y = 2 and so on. 
Thus there are q v nodes at level y, each with one connec- 
tion to a node at level y—1 and q nodes at level y + 1. The 
nodes at each level are partially connected to each other. 
Suppose that each node is connected to p other nodes at 
the same level. Note that p can be smaller than unity, 
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since nodes can represent groups of people, rather than 
individuals. 

As previously mentioned, each node x has a state s x , 
taking one of the values +1 or —1. We assume that the 
state is the same for all nodes at the same level. This is a 
mean field assumption which replaces individual states by 
the average state of their classes (intuitively, this comes 
from the structural symmetry of the nodes within each 
class. Sociologically, this is the same as the behavioral 
convergence of class members in Marxian theory, as men- 
tioned in the introduction section. Also note that, with 
this assumption, the tree graph collapses into a one di- 
mensional chain, with each node representing a class). 

We denote the strength of the links connecting level y 
to level y + 1 by J{y —> y + 1}, which is taken to be 
positive (the negative case implies anti-conformity which 
is not considered here, but might be an interesting line of 
study for future work). Define 



J{y-i^y} = 
J{y->-y} 

J{y^y+l} _ 
J{y^y} 



P 



(1) 



The flipping probability of each node is proportional to 
the weighted fraction of its neighbors with the opposite 
state. Denote the flipping probability for nodes at level 
y by w y . This means that with probability w v , we have 
s y (t + At) = —s y (t), and with probability 1 — w y , we have 
s y (t + At) = s y (t). Denote the set of nodes adjacent to 
the node at level y by N y . We have: 



^-JX 



J{x -t y}s x 



(2) 



which simplifies to: 



U'u 



qa + lxf3+px 1 



(fiSy-l + qaSy + 1 +pSy) 



(3) 

Hereafter, s y (t) represents the expected value of the state 
at level y at time t. Using the flipping probability, we find 



the time evolution of s y (t): 

qa + 13 
s„ =- I 7T-. — I s 



qa + j3 + p 1 " 

P 



qa + (3 +p 



Sy-l + 



qa 



qa + f3 + p 



)*v+i- ( 4 ) 



A. Continuous Approximation, Semi-infinite 
Structure 

First consider the case that the structure is limited 
from above but unlimited from below, i.e., y € Z + = 
{0,1,2,...}. A continuous approximation of Q can be 
obtained as a differential equation in y, by relaxing y to 
take values in M+ = {x G M : x > 0}. Note that || can 
be written equivalently as 



qa 



" \qa + 0+p 
qa- (3 



s qa + j3 +p 
In the continuous case we get 
qa 



(Sy + 1 + Sy-X - 2Sy) (5) 

(Sy-Sy-l)- (6) 



d t Sy(t) 



d L s v(t)+ ( ^ a T Z ) 5 9 s sW- 



(7) 



qa + P+pJ yy vV ' \qa + P+p 



This turns into the standard heat equation [M] for cf>(y, t) 
using the substitution: 



/(*) = <^(j/,*)exp 



— / n — \t y 

4qa(qa + p + p) 2qa 



(8) 

If initial state of the nodes are s y (0) and the head node is 
at some arbitrary state So(t), and defining A = qa + j3 +p 
for brevity, we get 




Sy(t) 



where we have defined 



^(y,t) + ^(y,t), (9) 



Hy,t) 



del' 



def 



(t = 0) exp 




exp 



(aq - (3) 



exp 



(y -0 s 

Atqa 



(A) 



exp 



ttqa 



(A) 



f3 — qa 
2qa 



-y 



s (t-X) 



exp 



(qa-(3) 2 (\) y 2 (A) 



Aqa(A) 



4qaX 



d\. 



(10) 



The first term in (|9| is transient and is due to initial con- 
ditions. The second integral can be approximated asymp- 



totically for long times: 



Sy(t) 



1 \qa-P\ 
qa + f3 +p 



exp 



( qa-P \ 
V 2ga J 



so[t-y 



qa + f3 +p 
\qa-P\ 
(11) 
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The condition on t ensures that the integrand has a max- 
imize^ and the asymptotic approximation is valid. In- 
tuitively, it implies that the "wave" of influence hasn't 
reached y at that time, and thus the expected state is still 
zero, assuming that function So(t) is turned on at t = 0. 

For the problem at hand, the head node is fixed at state 
so(t) = +1, for all t > 0. This implies that the function 
so(t) is the Heaviside step function u(t), resulting in 



exp 



( qa-p \ 
\ 2<?Q J 



Vv 



u t 



V- 



qa + f3 + p 
' \qa-P\ 



(12) 



Depending on the sign of ( qa a 13 ^ , the exponential factor 

can either grow or decay in y, giving rise to two different 
behaviors (Figures[2]and|3|. So there are two distinct cases 
to distinguish: 

• qa < j3 — > eventuates in s = +1 for all nodes. 

• qa > f3 — > the expected value of the spin decays in 



y, the grip of the head on the system debilitates. 



B. Continuous Approximation: Finite Structure 

Now we return to equation Q for the case of a finite 
hierarchy, y max — L. Let us define: 



a 



qa 
qa+/3+p 

qa+/3+p 



Using these, we rewrite Q: 



d t s y (t) = ady y Sy(t) + bd y Sy(t). 



(13) 



(14) 



The head node is fixed at So(t) = +1. We turn our 
attention to the case where there is complete disagreement 
between the head and the universal class. This means that 
the bottom- most class is fixed at si,(t) = — 1. 

With the substitution s y (t) — </>(y, t) exp {^y — 

turns into the standard homogeneous heat equation 
for (f>(y,t). Thus we arrive at: 



J 



2,t ^ _ ,nxy S 1-exp 



n sin I 



V L 



(111 .2 \ 

7r n a I b \ 
L I" 4a J 



4a 



exp 



(15) 



as depicted in Figures [4] and [5] For the steady state, when the time evolution vanishes, we have 



lims y (t) = a-^Y^ 



n=l 



x ^P (-g) +(-!)" exp (^l) 



4a 



(16) 



Observe that the steady state can be expressed in a sim- 
pler form by noting that ( 14 ) reduces to 



d 2 d 
a-r^Sy + o—Sy = 0, 
dy 2 y dy " 



(17) 



whose solution in the case of sx(t) = —1, s (t) = 1 is 

, , 2exp(^y) - 1 -exp(^L) t 
t-xx> i - exp {-f-L) 

The results are illustrated in Figures [6] and [7] 



Note that, in the case of agreement between the two 
ends, we have 



s y = +1 Vy, 



(19) 



which signifies the importance of the role of the bottom 
class. 



III. SECOND MODEL: CONSTANT FRACTION 
OF INTRA-CLASS LINKS 

In the second model, the number of links between the 
nodes of the same class is a constant portion of its popula- 
tion. This means that p y — rq v for some < r < 1. This 
is the case if we take the connection probability of pairs of 
classmates to be the same for every class. Generally, this 
case develops intuition about any structure in which ev- 
ery class can be modelled by the same graph topology, for 
which the average intraclass degree can be approximately 
taken to be linear in the number of nodes in the class. 

The evolution of the expected values of the states be- 
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time 



FIG. 2: Time evolution of the state on different sites on a semi-infinite structure when f3 = 1.2aq > ctq (thus the head node takes over), for 
y = 2 (top one) to y = 7. It is clear that, the closer to the head the class is located, the faster it complies. We simulated a chain of length 
5000 to approximate the semi-infinite chain for the y < 250 region, and averaged over 10000 Monte Carlo simulations. 
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FIG. 3: Steady state values on a semi-infinite structure when (3 = 0.7ctq < ctq, thus the influence of the head node dies out. We simulated a 
chain of length 5000 to approximate the semi-infinite chain for the y < 250 region. Monte Carlo simulations were added until the absolute 
value of the relative change in the mean state, averaged over all cites, became less than 10 -5 . 



qa + /3 + rqy) Sy [qa + (3 + rqV } Sy 1 



qa 

+ { qa + P + rqv) Sy+1 ' 



(20) 



Let us first assume the case of strong hierarchy, i.e., 
qa <C /3,p y ■ For a = 0, there is complete obedience to 
the superior class, and no weight is assigned to the inferior 
classes. For this case we have 



d t s y (t) 



+ rqV 



dySy(t). 



(21) 



If the head node is fixed at s = +1, which means s y —o(t) = 
u(t) then we have 



s y (t) = u 



t-y- 



q y - i 

Plnq 



(22) 



which implies that for long times eventually all levels con- 



form to the head node. 

The other polar case would be where qa 3> /3,p y which 
gives 



d t Sy(t) = dl y Sy{t) + dySy{t), 

whose solution for long times follows 

P -y 



Sy(t) 



t(*-y), 



(23) 



(24) 



which decays in y. 



For the steady state, ( 14 ) can be written as: 

\ d 2 s ( qa- /3 \ = Q 

qa + (3 + rqV J dy 2 Sy \qa + (3 + rqV ) dy Sy 



whose solution is of the form 

'(3 - qa 



K\ exp 



qa 



-II 



(25) 
(26) 
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FIG. 4 : Time evolution of states on a finite chain when f3 = 0.7ctq < ctq, implying the dominance of the universal class. The solid lines 
represent the theoretical prediction. The height of the structure is L = 40 in this simulation. The two at the bottom (y = 30, 35 ) are frozen 
at a state very close to — 1, because of the proximity to the universal class at y = L = 40. The results are averaged over 12000 Monte Carlo 

simulations. 
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FIG. 5: Time evolution of states on a finite chain when f) = l.2aq > aq, whose ramification is the dominance of the head. The solid lines 
represent the theoretical prediction. The height of the structure is L = 40. The one at the bottom corresponds to y = 35, at the vicinity of 
the universal class, but frozen near zero rather than —1. The state of the top two (y = 5, 10J, neighboring the head, is almost exactly +1. in 

this simulation. The results are averaged over 12000 Monte Carlo simulations. 



where K\ and K% are constants and are uniquely deter- 
mined for given boundary conditions. 

The solution in the case of s(L) = —1, s(0) = 1 coincides 
with ( 18 ), which does not depend on the value of r. Also, if 
s L (t) 



so(t) — +1 for all t, then we have lim 4 _ 
for all y. Note that, 



Sy(t) = 1 



• The value of r does not affect the steady state, but 
only how fast the system reaches it. 



• The only determinant factor of conformity is the 
relative value of j3 and qa, if the latter is greater, 
the state imposed by the head node decays, and it 
takes over otherwise. 



IV. GENERAL STEADY-STATE DISCRETE 
SOLUTION: FINITE HIERARCHY 



In this section we consider the general case for intraclass 
connectivity. We also take y to be discrete. Let the average 
degree of the class at level y be p(y). Then the evolution 
of the expected value of the state is 



Sy(t) = - 



qa + /3 


qa + /3 + 


p(y) 


P 




(?a + /3 + 


p(y) 


qa 




qa + P-r 


p(y) 



Sy(t) 



Sy-l{t) 



S y+1 (t). 



(27) 



Now we focus on the steady state, thus we drop the time 
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FIG. 6: Steady state on a finite chain with length L = 40 when f) = 0.7ctq < ctq . Note that the majority is converted to the lead of the 
bottom class. The results are averaged over 15000 Monte Carlo simulations, when the absolute value of the relative change in the mean 

state, averaged over all cites, became less than 10~ 5 . 




FIG. 7: Steady state on a finite chain with length L = 40 when fj = 1.2ctq > aq. The head node prevails clearly. The results are averaged 
over 8000 Monte Carlo simulations, when the absolute value of the relative change in the mean state, averaged over all cites, became less 

than 10 -5 . 



arguments. The steady state satisfies 





qa + 
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qa 




-p{y) 








qa 




-p{y) 




qa 




qa 




-p(y) 



Sy-l 

s y+1 = 0, (28) 

\</« t h t p\y> / 
or equivalcntly, 

- (aq + p)s y + fisy-x + aqs y+ i = 0, (29) 

which is independent of p(y). We emphasize that, the form 
ofp(y) does not affect the steady state, but merely how fast 
the system gets there. . 



Equation (29) is a linear difference equation with con- 
stant coefficients [65] ■ Therefore its solution has the form 
s y = Ki\\+K2\\. Coefficients K\ and K-i are determined 
using the boundary conditions. Ai and A2 are readily de- 
termined by plugging the ansatz s y = X y into ( 29 ) , which 



gives: 



aqX 2 - {aq + j3)X + [3 = 0. 



(30) 



The two solutions to this quadratic equation yield Ai , A2 . 



But to solve equation ( 29 1 , we require two boundary condi- 
tions to extract K\ and K2 ■ Thus we assume three distinct 
cases. 



• s(head) = +1 and s (bottom) 
the steady state we have: 



T. In this case for 



(JLY-2UY 

(««) 1 



(31) 



which is a transition between the two states. Major- 
ity is held by the head or the bottom if /3 is bigger or 
smaller than aq respectively. Figures [5] and [9] are de- 
pictions of these two possibilities for the case where 
p(y) varies exponentially with y, as mentioned previ- 
ously, meaning that a constant fraction of neighbors 



8 



from each class are linked. 



V. SUMMARY AND DISCUSSION 



s(head) = +1 and the bottom class is dangling, i.e., 



it follows (271, implying imitation of its upper class. 



In this case, we have 



aq / 



(32) 



Note that this cannot exceed unity. Divergence 
means that s = +1 is moving downwards. So we 
observe that: 



j3 > aq — > s(y) = 1 Vy. 



(33) 



Figures 11 and 10 illustrate this case. 



s(head) = +1 and s(bottom) = +1. In this case, 
limt-j.oo s(y) = 1 for all y, irrespective of values of 
/?, a and q. 



The implication of the latter is straightforward. The head 
node does not require to exert influence on the middle 
classes, but merely getting the universal class aligned suf- 
fices to dominate the entire hierarchy. This signifies the 
importance of media, or any other means by which this 
alignment could be managed. 



We studied the dynamics of binary states, for nodes on 
a pyramid-shaped hierarchical graph. The head node was 
given a fixed state, trying to impose it on the whole net- 
work. The conditions that facilitate or hinder the domi- 
nance of the head node were found. If classes assign more 
weight to classes below them, dominance of the head node 
is hampered. It was observed that in the case of discord 
between the head node and the bottom class in a finite 
hierarchy, one of them wins the majority. Another fact 
was that, the density of the intra-class links does not af- 
fect the steady state, but merely determines how quickly 
the system reaches it. Also, if the head node and the 
bottom-most class are fixed at the same state, the whole 
hierarchy will eventually comply, under any conditions for 
the inter/intra-class links. 

Possible extensions of this work include using a more 
realistic topology for the underlying graph to model the 
social structure, considering continuous and/or multidi- 
mensional states, deploying more than one classification 
attribute (e.g. Weberian class theory [57]), and having 
level-dependent ratios for weights assigned to lower and 
higher classes. Also, one can examine the effects of adding 
sparse long-range interactions to the hierarchy. Another 
modification would be adding fractions of stubborn agents 
within the middle classes and studying their effect on the 
final state. Also, one can add exogenous influence to the 
whole system, in the form of an external field, with differ- 
ent exposure degrees for different classes. This can model 
the effect of mass media. 
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